Effects of Coulomb interaction and tunneling on electron transport 
in coupled one-dimensional systems: from ballistic to diffusive regime. 
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A linear theory of electron transport is developed for a system of two ideal quantum wires, of length 
L, coupled by tunneling and Coulomb interaction. The interaction of electrons with acoustical 
phonons is included and the results are valid in both the ballistic and diffusive regime. In the 
ballistic regime, both tunneling and Coulomb drag lead to a negative transresistance Rtr, while in 
the diffusive regime the tunneling opposes the drag and leads to a positive Rtr- If L is smaller than 
the phase-breaking length, the tunneling leads to interference oscillations of the resistance that are 
damped exponentially with L. 
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In submicromcter-long quantum wires and at low tem- 
peratures the electron transport is mainly ballistic 1 and 
the wire conductance reaches its fundamental value of 
Go = e 2 /irh. In contrast, for sufficiently long wires 
this transport is limited by scattering processes. The 
tunneling of electrons between parallel quantum wires 
and/or interlayer electron-electron (e-e) interaction are 
essential as demonstrated, e.g., by experimental and the- 
oretical works on electron transport 2-13 along the layers 
in coupled double- wire systems. These works are mostly 
devoted to studying interlayer tunneling in the purely 
ballistic regime 2-10 or momentum transfer between the 
wires (Coulomb drag) 11-13 . The drag effect has been 
studied both in the diffusive 11 and ballistic 12 transport 
regimes as well as when the electrons are described by a 
Luttinger liquid 13 . 

Despite this progress, the description of electron trans- 
port in coupled quantum wires is substantially insuffi- 
cient. Even within the concept that the electrons are 
described by a normal Fermi liquid, two important ques- 
tions arise. The first one is how to describe the trans- 
port properties when both tunneling and the interac- 
tions of electrons with each other, and with impurities 
or phonons, are essential. The second question is how to 
bridge the gap between the ballistic and diffusive regimes 
in such a description. In this Letter we present a linear- 
response theory of electron transport in coupled quantum 
wires that gives a reasonable answer to both questions 
and the new results mentioned in the abstract. 

From kinetic theory to local description. We consider 
two parallel, tunnel- coupled ID layers of degenerate elec- 
trons adiabatically contacted to equilibrium reservoirs 
(leads) as shown and labelled in Fig. 1. We start 
from the quantum kinetic equation dp/dt + (i/h)[H + 
He + H e - p h,p] = for the density matrix p and as- 
sume that electrons interact with the Coulomb field (He) 
and with acoustical phonons(77 e _ p ^). Elastic scattering 
is neglected, i.e., we assume ideal wires. We employ a 
tight-binding description with basis that of the isolated 



left (I) and right (r) layer states in which the potential 
energy matrix is h = (A/2)<7 2 + Ta x when only the low- 
est subband is occupied in cither layer. Here &i are the 
Pauli matrices, A is the level splitting energy, and T the 
tunneling matrix element. 

The kinetic equation can be written 14 as one for the 
Keldysh's Green's function G h . Below we consider the 
case when the characteristic spatial scale of the electronic 
distribution is large in comparison with the Fermi wave- 
length wh/pF and use the Keldysh's matrix Green's func- 
tion in the Wigner representation G~ + ^(p, x), where p 

and e are the momentum and energy. In the steady-state 
regime the time-averaged Green's functions are linearized 
in the manner Gf(p, x) = G ( ° )ap {p) + SGf{p, x), where 
a and stand for + or — . The unperturbed part 
Gg °^(p) is expressed through the retarded and ad- 
vanced matrix Green's functions and the equilibrium 
Fermi-Dirac function /(e) = 1/[1 + e (e- ^/ fcBTe ] in the 
usual way 14 . The linearized kinetic equation reads 



~ {% ,^G- £ + (p, x) } + i [h, SGs + (p, *)]+#, Gf + (p)\ 
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~(p)} = iSI(e,p, x). 
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Here {...} denotes anticommutators, v p = PiVi p + P r v rp 
is the^ diagonal matrix of the group velocities, and Pi 
and P r are the projection matrices. We consider only 
the case of equal group velocities in the layers, when 
vip = v rp = v p = p/m. Further, (p is the matrix of 
the self-consistent electrostatic potential resulting from 
the perturbation of the electron density. The collision 
integral is given as 14 X = -[E - +G++ + S -- G -+ + 
G -+ S++ + G -- S -+ ], where t a P is the self-energy due 
to Coulomb and electron-phonon interactions, and where 
the arguments of all functions are e, p, and x. This cor- 
responds to a quasiclassical description of the scattering. 
Despite the approximations made, Eq. (I), due to its 
matrix structure, is not reduced to a classical Boltzmann 
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equation. 

We assume that the Fermi energy is large in compar- 
ison with both T and A. We sum up Eq. (1) over the 
electron momentum p in the regions of forward (+) and 
backward (— ) group velocities and introduce the non- 
equilibrium part of the energy distribution function ge{x) 

gf(x) = (l/27ri) Jjp\v p \5G^ + (j>,x). (2) 

Since SG ^ is essentially nonzero only in narrow intervals 
of energy and momentum near the chemical potential fi 
and Fermi momentum, we can replace \v p \ by the Fermi 
velocity vp common to both layers. The result is 

±v F dgi(x)/dx + (i/h)[h,gf(x)] = 51±{e,x) (3) 

where SI±(e,x) = (l/2nh) J ± \v p \dp ST(e,p,x) depends 
on both g + and g~ since it accounts for both forward- 
and backscattering processes. Notice that the absence of 
(p in Eq. (3) is not a simplification or an approximation: 
it results exactly from the integration of Eq. (1) over 
momentum. 

The matrix kinetic equation (3) is equivalent to eight 
scalar equations for the four components of g + and those 
(four) of g~ . The boundary conditions for them arc 
determined, in the Landaucr-Biittiker-Imry approach, 
by the distribution functions of the leads described 
by four chemical potentials /in, /ii r , [i2i , and [i2 r , 
so that <7e"(0) = — (df(e)/de)[Pi5fj,u + P r 8ni r ], and 

g~e{L) = —{df(e)/de)[PiS^2i+PrS^2r]- The forward- and 
backward-propagating states are "connected" to leads 
"1" and "2", respectively. The nondiagonal components 
vanish at the contacts since the tunneling is absent be- 
yond the region x = [0, L\. 

We now introduce the local chemical potential matrix 
5fi k {x) = PiS^f(x) + P r 5nt(x) + a x 6fj£{x) + a y S^f(x) 
according to 5fi ± (x) = j de gf(x). Assuming that 
backscattering occurs much less frequently than forward- 
scattering, which is the case when the e-e scattering dom- 
inates over the electron-phonon one, and that tunneling 
also occurs much less frequently than forward-scattering, 
we can write [ffe{x)\u,rr — ~(9f{£)/de)S/j^ z r (x). Below 
we omit the symbol "5" in 6jl (x) and in the potentials 
of the leads since all potentials are measured from the 
same equilibrium value \i. Then integrating Eq. (3) over 
the energy we obtain eight coupled, first-order differential 
equations: 

± dfif /dx + {pf - fif )(l/l P + l/l D ) 

-(nt-fi)/l D -2t F tf=0, (4) 

±dfif/dx + (pf - $){1/Ip + l/l D ) 

-(fif -fif)/l D + 2t F ^ =0, (5) 

±d^/dx + S F Hy + nt/h = 0, (6) 
±dy±/dx - 5 F nt + Vy/te + t F {nf — fif) = 0. (7) 



Here tp = T/Tivp,5p = A/hvp. The evaluation of SI 
was carried out assuming a weak tunnel coupling, when 
T is small in comparison to the imaginary part of the self- 
energies. The characteristic lengths lp 7 Id and lc, result- 
ing from the collision integral 5I±(s, x), are expressed, 
respectively, through the phonon-assisted ID transport 
time 15 rp, the ID Coulomb drag time 11 td, and the 
phase-breaking time tc (which describes the suppression 
of tunnel coherence) as lp — 2vpTp, Id — 2upTp>, and 
lc = vpTc- All these times are microscopically justified. 
For tc, if we account only for the lowest-order Coulomb 
contributions, we obtain the new result 

1/t c [e 4 ln 2 (w/a)A/2irh 3 e 2 vp]cotanh(A/4:k B T e ) (8) 

where e is the dielectric constant, w is the distance be- 
tween the wires and a the wire width. Estimating tc 
from Eq. (8), we find lc "C lp and lc <C Id, because of 
the weakness of the electron-phonon coupling and of the 
e-e backscattering, respectively. 

Equations (4)- (7) with the boundary conditions 
/"/"(O) = Miii A*r( L ) = A*2Z, M+(0) = fj, lr , li~(L) = [l 2r 
and nt,y{fy — Vx,y(L) — give a complete description 
of electron transport in double quantum wire systems in 
a wide range of regimes starting from the purely bal- 
listic regime, for L <^ lc, to the diffusive regime, for 
L ^> Ip,Id- The local currents flowing in the layers 
j = l,r are expressed by the Landauer-like formula 

J j (x)=G [n+(x)-vj(x)]/e. (9) 

Below, to characterize the effects of drag and tunneling, 
we calculate the transresistance Rtr = [fJ-u — A 1 2i]/eJ r , 
where J r is the current injected in the wire r when no 
current is allowed to flow into wire I. This is the typical 
setup for the drag experiments. In this case J r (0) = 
J r (L) = J r and J;(0) = Ji{L) = 0. We also calculate the 
direct resistance R — [fi\ r — H2r]/eJ r - 

Results. The solution of Eqs. (4)-(7) is obtained for 
all potentials in the form of A + Bx + ^2±Cfe ±XiX , 
where Af are the roots of a cubic equation. Explicitly, 
for lc <C Ip,Id and weak tunnel coupling tp <C Iq 1 , 
the three roots are Ai = A ~ 2(1/1^ + l/lpli) 1 ^ 2 and 
A 2 , 3 = A± ~ l/l c ± iS F , where 1/h = l/l P + 2/l D 
and lp = vpTp- The tunneling time Tp is given by 
1/ Tt = (1/t c )2T 2 /(A 2 + H 2 /t c ). The terms with A 
describe long-scale variations of the chemical potentials, 
while those with A± correspond to short-scale variations. 
Two cases can be considered. 

Long wires, L ^> lc ■ The solutions containing A± arc 
not essential in the calculation of the currents. Consid- 
ering only the solutions with A, we find 

R = (7rft/2e 2 )[2 + L/l P + y/l + l T /h tanh(AL/2)], (10) 
R TR = (nh/2e 2 )[L/lp - y/l + l T /h tanh(AL/2)]. (11) 

In the ballistic regime [L <^ lp,lr) we obtain the usual 
resistance R~ nh/e 2 — Gq 1 , while the transresistance is 
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given by Rtr — — {~Kh/e 2 )L[l/lD + l/2lT\- As seen, Rtr 
is small, always negative, and proportional to L. If we ne- 
glect tunneling, the resulting Rtr describes the Coulomb 
drag in the ballistic regime 12 . In the diffusive regime, 
i.e., when L increases and the relation L ^> lp holds, 
we obtain R ~ (nh/e 2 )L[l/l P + l/l D ] for XL/2 <C 1. 
This resistance, if one omits the drag contribution, is ex- 
pressed in terms of the Drude conductivity a = L/ R = 
e lp/irh = e 2 nTp/m, where n is the ID electron density 
in layer r. Then R TR = -~(irh/e 2 )(L/l D )[l - (L/L ) 2 } 
where L$ = (<o1 2 p It /Id) 1 ^ 2 ■ Introducing the drag transre- 
sistivity it%/e 2 lp> and the ID-ID tunneling conductance 
e 2 piD/vr = 2e 2 /-kKIt, where pw is the ID density of 
states at the Fermi level, we formally obtain the result of 
Rcf. 16 where a competition of drag and tunneling was 
investigated, for double quantum well systems, only in 
the diffusive regime. As the factor [1 — (L/L ) 2 ] shows, 
in this regime the tunneling opposes the drag and the 
transrcsistance increases, changing its sign from negative 
to positive. If XL/2 ~ 1, Rtr is large and always compa- 
rable to R, because a considerable fraction of the current 
penetrates in the I layer due to tunneling. In Fig. 2 we 
show the length dependence of Rtr given by Eq. (11) 
for different relative contributions of the Coulomb drag 
and tunneling. Rtr is negative for small L but it always 
changes sign and becomes positive as L increases. This 
occurs at smaller L when the tunneling is stronger (larger 
Ip/It) or the drag is weaker (smaller Ip/Id)- 

A peculiar transport regime, without backscatter- 
ing, can be realized in tunnel-coupled magnetic edge 
states 17 ' 18 . Assuming l/l P = l/l D = in Eqs. (10) 
and (11), we obtain the result of Rcf. 18 in the 
form R = (nh/e 2 ) [1 + (1/2) tanh(L/Z T )] and R TR = 
-(7r?i/2e 2 )tanh(L/; T ). 

Short wires L ~ lc <C Ip,Id- Electrons pass through 
the wires almost without backscattering, R ~ nh/e 2 , and 
Rtr is small. However, an electron tunneling between 
the layers does not lose its phase memory completely; as 
a result tunnel coherence effects can take place and give 
additional contributions to R and Rtr- 

where 

$(L) = (l c 2 + 4)- 1 [2((5 F //c)e" L/ ' c sm(5 F L) 

+ (l C 2 - - z~ L/lc cos(S F L))}. (14) 

The term proportional to &(L) describes oscillations 
of the resistance damped due to the factor exp(— L/Iq)- 
Its relative contribution to Rtr is not small at L ~ lc- 
The periodic behavior can be described as a result of 
tunncling-assisted interference of electron waves of the 
left and right layers on the coupling length L. Due to 



a finite A, these waves have slightly different phase ve- 
locities. Figure 3 shows the dependence of $(L) as a 
function of SfL; we used Eq. (8) at zero temperature so 
that 1/lc = k\Sf\, where k is a dimensionless parameter. 

The level splitting A can be changed not only by a 
transverse voltage, but also by a magnetic field 9 H per- 
pendicular to the plane of the quantum wires. For suffi- 
ciently weak H the results presented so far still hold with 
the phase SfL in Eq. (14) changed to S P L + 2ir4>/(j)o, 
where <po = h/e is the magnetic flux quantum and 
4> = HwL the flux enclosed by the area between the 
wires. Though the double-wire system does not form a 
closed current loop, this should lead to Aharonov-Bohm- 
type oscillations in R and Rtr- 

In very short wires, with L <C IcSp 1 , we have Rtr — 
{-Kh/e 2 )[—L/lo — L 2 t 2 F /2\: in contrast with the drag con- 
tribution, that of tunneling shows a L? dependence. 

Finally, we apply Eq. (3) to the purely ballistic regime 
in which we can neglect the collision integral and need not 
make the assumption of weak tunnel coupling. Then we 
obtain Eqs. (4)-(7) without the terms containing lp, lp, 
and lc', they describe oscillations of the electronic wave 
packets between the layers due to coherent tunneling. 
Solving them, we obtain 

R = (nh/e 2 )[l - rsin 2 ?A][l - 2r sin 2 t/-]" 1 , (15) 
Rtr = -(nh/e 2 )rsm 2 ^[1 - 2r sin 2 V-]" 1 , (16) 

where r = 2T 2 /A 2 T , A T - (A 2 + 4T 2 ) 1 / 2 and V = 
ApL/2hvF- The periodicity of the oscillations becomes 
the same as in Eqs. (12) and (13) if one replaces A by At- 
However, since the tunnel coupling is not weak, the oscil- 
lations described by Eqs. (15) and (16) have large ampli- 
tudes. In particular, when A is small, both R and Rtr 
show giant oscillations with amplitude large in compari- 
son with 1 /Go- The experimentally observed 6 resistance 
oscillations, resulting from tunnel coupling over a finite- 
length region, in ballistic quantum wires were of small 
(~ 0.5 Kf2) amplitude. This is not surprising because 
there are many factors, e.g., long-scale inhomogeneities, 
inelastic and elastic scattering, which compete against 
the tunnel coherence. Previous theoretical studies of the 
purely ballistic transport regime 2-5,7-9 were based on a 
quantum-mechanical calculation of the electronic wave 
transmission, while we recover the essential results from 
a quantum-kinetic analysis. 

In conclusion, we developed a linear-response, steady- 
state theory of electron transport in parallel onc- 
dimcnsional layers coupled by tunneling and Coulomb 
interaction and contacted to quasi-equilibrium reservoirs. 
The quantum-kinetic description of the problem leads to 
linear differential equations describing the distribution of 
local chemical potentials of the forward- and backward- 
moving electrons. The results for the direct resistance R 
and the transrcsistance Rtr are valid in the whole range 
from the pure ballistic to the diffusive regime. Rtr al- 
ways depends nonmonotonically on the wire length L: in 
the ballistic regime tunneling and drag lead to a negative 
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transresistance Rtr, while in the diffusive regime the 
tunneling opposes the drag and leads to a positive Rtr- 
When L is smaller than the phase-breaking length, the 
tunneling leads to interference oscillations of the resis- 
tance damped exponentially with L. 
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Fig.l O.E.Raichev and P.Vasilopoulos, Effects of... 




Fig. 2 O.E.Raichev and P.Vasilopoulos, Effects of. 



